Analysing J/^ Production in Various RHIC Interactions with a 
Version of Sequential Chain Model (SCM) 



in 



0^ . p. Guptaroy^* Tarun K. Garain^, Goutam Sau^j 

^ ! S. K. Biswas^f & S. Bhattacharyya^^ 

\ ^ Department of Physics, Raghunathpur College, 

513 ' Raghunathpur 723133, Purulia, West Bengal, India. 

, ^ ■ ^ Department of Mathematics, Raghunathpur College, 

Raghunathpur 723133, Purulia, West Bengal, India. 
^ Beramara Ram Chandrapur High School, 
South 24-Pgs,743609(WB),India. 
f~| ■ West Kodalia Adarsha Siksha Sadan, New Barrackpore, 

Kolkata-700131, India. 

flLi; ^ Physics and Apphed Mathematics Unit (PAMU), 

^ ■ Indian Statistical Institute, Kolkata - 700108, India. 



> 

(N 
(N 

\^ • We have attempted to develop here tentatively a model for J/^* production in p+p, d+Au, 

\ Cu + Cu and Au + Au collisions at RHIC energies on the basic ansatz that the results of 

Q*^ ■ nucleus-nucleus collisions could be arrived at from the nucleon-nucleon [p + p)-interactions 

with induction of some additional specific features of high energy nuclear collisions. Based 
on the proposed new and somewhat unfamiliar model, we have tried (i) to capture the 
properties of invariant p^-spectra for J /"^ meson production; (ii) to study the nature of 
centrality dependence of the p^-spectra; (iii) to understand the rapidity distributions; (iv) 
' to obtain the characteristics of the average transverse momentum < px > and the values of 

< p^ > as well and (v) to trace the nature of nuclear modification factor. The alternative 
approach adopted here describes the data-sets on the above-mentioned various observables 
in a fairly satisfactory manner. And, finally, the nature of J/^t-production at Large Hadron 
Collider(LHC)-energies deduced on the basis of our chosen model has been presented in a 
predictive way against the RHIC-yields, both calculated for the most central collisions and 
on the same model. 
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"Finally J/^' suppression, which for more than 20 years has represented the gold-plated signature 
of deconfinement, is not understood" - M. J. Tannenbaum [l]. 

1 Introduction and Background 

At the very start let us first address two somewhat imposed but important and interesting 
questions: (i) why are we so seriously interested in persuing the ^'-studies at both RHIC and 
LHC?; (ii) Why do we turn to a somewhat new model, leaving aside the several versions of the so- 
called standard model? In fact, the much- valued observation contained in and captivated by the 
quotation at the top from a work by Tannenbaum [1] answers, in part, simultaneously both the 
questions. But, an elaboration on the background and the perspective, we feel, is quite necessary 
for better understanding and comprehension of both the points to the desired degree and depth. 
Why has the J/^'-suppression attracted the most attention as the likely "gold-plated" signal 
[2]? The pioneering work of Matsui and Satzp] established at the very beginning a few points: 
(i) as a hard QCD process, the heavy charm pair production takes place very easily, (ii) the 
Debye screening of the QGP prevents formation of a J/^' state in heavy ion collisions, (iii) the 
low temperatures at the hadronization do not really permit of the proposed charm anticharm 
pair kinematically; and quite soon after them, it was further proposed that the suppression 
pattern ought to have a characteristic transverse momentum dependence. Later, perturbative 
calculations based on perturbative Quantum Chromodynamics (pQCD) established the fact that 
both the suppression signal itself as well as its pT-dependence could be mimicked by the mundane 
nuclear shadowing. It was thus quite evident from the early days that a detailed quantitative 
analysis would be necessary to disentangle the effects of Debye screening in QGP. It has since 
been recognized that other effects, notably the simple absorption of the produced J/^* in the 
nucleus, cause suppression of the produced J /'^ in all nucleus-involved collisions, categorized 
physically as normal suppression. The additional or so-called 'anomalous' suppression as 
the possible signal of QGP, shall have to be investigated into only after accommodating the 
normal suppression. But, the estimation of these proposed normal/anomalous suppressions 
got absolutely entangled in multiple problems in pQCD due to lack of reliable and consistent 
knowledge about the several related factors. The reasons are explained quite well by Bhalerao 
and Gavai [2]. Besides, a detailed review made by Rapp et al.[4j has pinpointed and elucidated 
the several conceptual, difficulties and experimental hurdles encounterde by the Standard Model 
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(SM) in understanding and interpreting the various aspects of the J/^'-data. So, in order to 
escape these problems and to minimise the difficulties due to uncertainties, we have tried here 
an alternative methodology for understanding the totality or a large of part of the data on 
J/^'-production in some ultrahigh energy collisions like p + p, d + Au, Cu + Cu and Au + Au 
interactions which comprise the totality of the contour of the experimental RHIC studies [5]- 
[7]. The present study aims at probing and comprehending some of the major properties and 
data-characteristics of J/^'-production in all the aforementioned RHIC-studies in a model-based 
manner. And the model of our choice is named to be the Sequential Chain Model (SCM) which 
has no QCD-tag and with some ancillary physical ideas. Very initially, and even now, the goal 
of high energy heavy ion physics has been to study 'quantum chromodynamics'(QCD) in a 
regime of high temperature, high density and large reactive volume. The hope was and is to 
find conclusive evidence that 'QCD' undergoes a phase transition at a critical temperature from 
a confined state, where so-called quarks and gluons are bound in colourless hadron states, to a 
hypothetical deconfined state named 'quark-gluon plasma' (QGP), where 'quarks' and 'gluons' 
can explore volumes larger than the typical hadron radius (r ~ Ifm). But in our endeavours 
made here we would try to avoid using all these standard phrases, propositions and methods; 
rather we would try to build up an entirely new and alternative approach to understand the 
production- mechanism for J/'^ production and interpret the relevant observables in the light of 
an alternative approach outlined in the next section. 

In one of our previous works [8], we demonstrated in no uncertain terms that the produc- 
tion of J/'^ particles was neither 'anomalously' suppressed, nor enhanced; rather it was "just 
normal" where the term 'normal' is to be interpreted as the similarity of the J/^'-behaviour 
with production characteristics of some other lighter mesons and hadrons. In the past quite 
consciously, we never differentiated between initial state and final state effects. Herein too we 
do the same. And finally now we will be probing here whether there could be any special effect 
arising out of a smaller system-size for Cu + Cu compared to that for Au + Au reactions. 

Our plan of work presented here is as follows. In the next section we try to present a model not 
based on the concepts of quarks and gluons; rather it is founded on altogether different facade 
presented in details in the Refs|9]-[T2j. In section 3 we give the results and a general discussion on 
them. And in the last section (section 4) we offer the final remarks which essentially summarise 
the conclusions of this work. 
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2 The Model and the Approach 



The description of the model-based features would be subdivided into two parts. The first part 
gives a brief overview of the J/^'-production mechanism in nucleon-nucleon (p + p) interaction 
in the context of the Sequential Chain Model (SCM). Thereafter, the relevant transition for 
different observables from p + p to A + B interactions will be discussed. 

The outline and characteristics of the model we use here are obtained, in the main, from 
one of our previous works [8]. According to this model, the photons released by two dominant 
varieties of vector mesons which are produced in a q-oj-tt sequential chain with q (rho) and 
cj (omega) alternating with pion(s). According to this model, called Sequential Chain Model 
(SCM), high energy hadronic interactions boil down, essentially, to the pion-pion interactions; 
as the protons are conceived in this model as p={i:~^ i:^ •&) ^ where is a spectator particle needed 
for the dynamical generation of quantum numbers of the nucleons [9]- [16]. The particle i? is 
called a spectator particle because it does not and cannot take part in the strong interactions, 
as it is identified to be a muonic neutrino and is taken to be a Majorana spinor. The internal 
quantum numbers like isospin, strangeness and baryon numbers can be related to the internal 
angular momentum of the constituents and we thus do have a geometrical origin of the basic 
SU(3) symmetry of the hadrons. And it is taken that these constituents move in a harmonic 
oscillator potential with orbital angular momentum 1/2/i in such a way that, the two values of 
the third component of the orbital momentum represent the two states of matter: particles and 
antiparticles. The multiple production of J/^-mesons in a high energy proton-proton collisions 
is described in the following way. The secondary vr-meson or the exchanged ^>-meson emit a 
free cj-meson and pi-meson; the pions so produced at high energies could liberate another pair 
of free q and trapped w-mesons (in the multiple production chain). These so-called free q and 
w-mesons decay quite a fast into photons and these photons decay into ^ or ^' particles, which, 
according to this alternative approach is a bound state of 0.0. or O'O' particles. The production 
of J/^'-mesons is shown schematically in Figure 1. The field theoretical calculations for the 
average multiplicity of the ^-secondaries and for the inclusive cross-sections of the psi-particles 
deliver some expressions which we would pick up from [8] and |17j . 

The inclusive cross-section of the \I'-meson produced in the p + p collisions given by 




exp( 



-5.35(p^ + m^/^) 





< nj/q, >l+p (1 - x) 
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where the expression for the average multiplicity for \I'-particles in p + p scattering is 

<nj/^>p+p = 4x10-^1/^ (2) 

In the above expression for eqn (1) the term |Cj/^| is a normalisation parameter and is as- 
sumed here to have a value = 0.09 for Intersecting Storage Ring(ISR)energy, and it is different 
for different energy region and for various collisions. The terms pT, x and rriji^, represent the 
transverse momentum, Feynman Scaling variable and the rest mass of the J/^' particle respec- 
tively. Moreover, by definition, x = Ipi^j \fs where pi, is the longitudinal momentum of the 
particle. The s in equation (2) is the square of the cm. energy. 

The second term in the right hand side of the equation (1), the constituent rearrangement 
term arises out of the partonic rearrangements inside the proton. But this proposed 'rearrange- 
ment' factor is different from what are implied by the terms 'parton recombination' or 'parton 
coalescence' in the standard literatures on ^'-studies. In fact, this is also different from what was 
meant by parton rearrangement by Bleibel et al. p8]. It is established that hadrons(baryons 
and mesons) are composed of few partons. At large transverse momenta in the high energy 
interaction processes the partons undergo some dissipation losses due to the impact and impulse 
of the projectile on the target and the parton inside them (both the projectile and the target) 
they suffer some forced shifts of their placements or configurations. These rearrangements mean 
undesirable loss of energy , in so far as the production mechanism is concerned. The choice 
of would depend on the following factors: (i) the specificities of the interacting projectile 
and target, (ii) the particularities of the secondaries emitted from a specific hadronic or nuclear 
interaction and (iii) the magnitudes of the momentum transfers and of a phase factor (with a 
maximum value of unity) in the rearrangement process in any collision. And this is a factor for 
which we shall have to parametrise alongwith some physics-based points indicated earlier. The 
parametrisation is to be done for two physical points, viz., the amount of momentum transfer 
and the contributions from a phase factor arising out of the rearrangement of the constituent 
partons. Collecting and combining all these, we proposed the relation to be given by [16] 

Nr = ^< Npart e, (3) 

where < Npart > denotes the average number of participating nucleons and 9 values are to 
be obtained phenomenologically from the fits to the data-points. In this context, the only 
additional physical information obtained from the observations made here is: with increase in 
the peripherality of the collisions the values of 9 gradually grow less and less, and vise versa. 
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is a phase term related to the geometry and number of the rearrangement of the partons of 
the nucleus. And the number of rearrangements is constrained by the impact parameter which 
decides the centrality/peripherality of the collisions reflected by Npart- 

The complex calculations of the Feynman diagram with the infinite momentum frame tech- 
niques give rise to such forms with some assumptions and approximations based on the bound- 
edness of px (transverse momentum)-values. For calculation purposes we took fgujn, d'yui and g^^ 
as the point-coupling of Fig.l. Moreover, in this context, we would have to reckon some other 
model-dependant features: (1) There are total nine probable diagrams for the various exchanges 
of vrvr — > vrvr. So the final contribution has a product term of 9^; (2) Secondly, the cm. energy 
of the TTTT system has to be translated into the actual p + p collisions for which one derives 
a conversion factor s,r7r = ^Spp. The previous paper [S] does not contain the product term 
(i.e. the factor) 1/p^^. This multiplier term has been introduced from physical considerations 
arising out of large-py (hard) reactions at RHIC energies, in fact, this term has been inducted 
absolutely phenomenologically. But the rest is an outcome of very rigorous calculations based on 
'soft' (small-p-r) production of particles. In nature 90% of the particles are produced with very 
small pt- And our initial focus during the period 1975-1980 was to concentrate on generalized 
particle-production characteristics in nature. 

In order to study a nuclear interaction of the type A + B Q + x, where A and B are 
projectile and target nucleus respectively, and Q is the detected particle which, in the present 
case, would be J/^'-mesons, the SCM has been adapted, on the basis of the suggested Wong 
|19| work to the Glauber techniques by using Wood-Saxon distributions |20j-|21j. The inclusive 
cross-sections for J/^ production in different nuclear interactions of the types A + B J/^ + X 
in the light of this modified Sequential Chain Model (SCM) can then be written in the following 
generalised form as [T3 ] -[l4 ] : 



E-t^Iab^j/^+x = Pj/^PT exp(-c(pf. + m]/^)) exp(-1.923< rij/^ >pp x). (4) 



the details of the calculations to be obtained from Refs.[13j-|14j. the set of relations to be used 



Here, in the above set of equations, the third factor gives a measure of the number of wounded 
nucleons i.e. of the probable number of participants, wherein Aas gives the probability cross- 




where Pj/qi, and c are the factors to be calculated under certain physical constraints. With 




1 + q(A1/3 + 51/3) 



1 



(5) 
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section of collision with 'i?' nucleus (target), had all the nucleons of A suffered collisions with B- 
target. And BaA has just the same physical meaning, with A and B replaced. Furthermore, cta 
is the nucleon(proton)-nucleus(A) interaction cross section, as is the inelastic nucleon(proton)- 
nucleus(B) reaction cross section and a^B is the inelastic AB cross section for the collision 
of nucleus A and nucleus B. The values of a^B, o'a, ctb are worked here out in a somewhat 
heuristic manner by the following formula [22] 

^Xl' = ^0 (<o,e.*.e + Al^r,et " (6) 

with fjo = 68.8 mb, 6 = 1.32. 

Besides, in expression (5), the fourth term is a physical factor related with energy degradation 
of the secondaries due to multiple collision effects. The parameter a occurring in eqn.(5) above 
is a measure of the fraction of the nucleons that suffer energy loss. The maximum value of a is 
unity, while all the nucleons suffer energy loss. This a parameter is usually to be chosen |19j . 
depending on the centrality of the collisions and the nature of the secondaries. 

The "-Pj/>ii" factor in the expression (5) accommodates a wide range of variation because of 
the existence of the large differences in the in the normalizations of the J/^ cross-sections for 
different types of interactions. The constraints for calculating -Pj/ij, are mentioned above. 

The parameter W/j' of Eqn. (4) has been calculated from the Eqn.(3). The constraints for 
calculating Nf(^ have been stated in the paragraph written below the Eqn. (4). 

The parameter 'c' is a energy and rapidity dependent term. It is clear from eqn. (1) that 
energy dependent factor comes from the average multiplicity of J/'f particles(~ s~'^'^^^)and the 
rapidity factor comes from (1 — x) where x = ^'"t smh _ g^^ ^ j^g^g (Jiffgj-gnt values for different 
energies and rapidity regions. 

The parameters used in this work are {'i)Cj/^ is a normalization term which is by nature 
species-dependent and energy-dependent; (ii) Nji = 4 < Npart >^^^ 6; this is determined by 
fits to the invarient cross-section data. It is naturally found to be both secondary-species and 
centrality-dependent; (iii) Finally, a is the fraction of the nucleus that suffer energy by under- 
going multiple collisions and this a-factor too has weak centrality dependence. 

The parameters are generally obtained from the minimization of of the experimental spec- 
trum for invariant cross-section which involves some rigorous statistical methods and complica- 
tions. But the main issues on parameter-constraining are the following few points: (i) availability 
of the large number of experimental data-points; (ii) smallness of the uncertainty-ranges in the 
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measured data. Unfortunately, for J/'^ production, in any high energy cohision whatsoever, 
the number of data points are relatively sparse, especially for the large transverse momentum 
values. Secondly, the uncertainty ranges in the measurements are also quite considerable, for 
which we completely dropped the idea and bid to devise the methods of constraining the evolved 
parameters. Very recently, Arleo and d'Enterria [23] studied on constraining the parameters for 
very simple and abundantly produced neutral pion production in pp-scattering at ^/s=22A GeV 
and ended up with a statement like " At high pxs the fit is completely unconstrained due to 
the lack of data, and its uncertainty is very large." 



3 Model-based Analyses and the Results 

Here, what we need to emphasize is that the observables are different at the same energy. So 
let us split up the cases and treat the data available on them, individually and on a case to case 
basis. 

3.1 J/^ (total) Crosssections \n p + p Interactions 

As the psi-productions are generically treated rightly as the resonance particles, the standard 
practice is to express the measured J/'^ (total) crosssections times branching ratio to muon 
or electrons, i.e. for lepton pairs , that is by Bu'ap!^ [23]. And this is usually done for all the 
symmetric nuclear collision like Au + Au, Ph + Pb and Cu + Cu etc. 

By using expression (4) we arrive at the expressions for the differential cross-sections for the 
production of J/^'-mesons in the mid and forward-rapidities (i.e. \y\ < 0.35 and 1.2 < \y\ < 2.2 
respectively) in p + p collisions at ^snn=200 GeV at RHIC. 

^i?//'^^lp+p-.j/M = 6.1p^^-^«3exp[-0.13(p| + 9.61)] for |y| < 0.35, (7) 

and 

1 ^ 



2^^^5n'^^|p+p-.j/M = 6.5p^i-i«^exp[-0.16(p^ + 9.61)] for 1.2 < \y\ < 2.2. (8) 

For deriving the expressions (7) and (8) we have used the relation x ~ ^p^™- = 2™,^ snihycm ^25\, 
where my , Vcm are the transverse mass of the produced particles and the rapidity distributions. 
mj/^ ~ 3096.9 it O.OllMeF [25] and Bn', the branching ratio is for muons or electrons i.e. its 
for lepton pairs J/^' fi'^ fi~ / e'^ e~ , is taken as 5.93 it 0.10 x 10"'^ [25] in calculating the above 
equations. 



Table 1: Values of ftj/vp and Npi for different centrality regions in Cu + Cu collisions 



Centrality [\y\ < 0.35] 


0-20% 


20-40% 


40-60% 


60-92% 




6.40 X 10-^ 


3.82 X 10-5 


1.22 X 10-5 


2.10 X 10"^ 


Nr 


0.829 


0.709 


0.625 


0.615 


Centrality [1.2 < \y\ < 2.2] 


0-20% 


20-40% 


40-60% 


60-92% 




6.38 X 10-5 


2.82 X 10-5 


0.82 X 10-5 


1.90 X 10-5 


Nr 


0.829 


0.709 


0.625 


0.615 



In Figure (2), we have drawn the solid lines depicting the model-based results with the help 
of above two equations (7) and (8) against the experimental measurements [26] .127]. 

3.2 Invariant yields of J/^P Particles in d + Au, Cu + Cu and Au + Au Collisions 

From the expression (4) , we arrive at the invariant yields for the J/ ^'-production in d + Au 
J/^ + X reactions for mid and forward-rapidities. 

^^^U+A«^j/M = 7.25 X 10-^p^°-''' exp[-0.13(p| + 9.61)] for \y\ < 0.35, (9) 
and 

1 d'^N 

^^^U+A«^j/M = 4.25 X lO-^p^O-629 exp[-0.16(p^ + 9.61)] for 1.2 < \y\ < 2.2. 

(10) 

In Figure (3), we have drawn the solid lines depicting the model-based results with the help of 
above two equations (9) and (10) against the experimental measurements |28j . 

For the case of Cm+Cu collisions at RHIC, the values of = Pj/,!, xexp(— 1.923< nj/^ >ppX^ 
and Nr of expression (4) are given in Table 1 for the rapidities |y| < 0.35 and 1.2 < \y\ < 2.2 
respectively and they are plotted in Figure 3. The exponential parts for different rapidities are 
remain the same as in equation (7) and equation (8). As previously stated, the exponents Nr 
of pt depend on the following three factors mentioned in the previous section. The values of 
< Npart > of equation (3) for different centralities of Cu + Cu collisions have been taken from 
|29j and the calculations have been done accordingly. The experimental results for the invariant 
yields of J/^ production as a function of transverse momenta [30] at different centrality values 



^The exponential part, here, contains s and x i.e. the squared cm. energy and the rapidity factor j/cm- This, 
in course of calculations, gives some fixed values and is multiplied with Pj/q, to give the final product bj/^. 
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Table 2: Values of hji^, and Nr for different centrality regions in Au + Au collisions 



Centrality [\y\ < 0.35] 


0-20% 


20-40% 


40-60% 


60-92% 




1.32 X 10^3 


0.86 X 10-3 


1.65 X 10"^ 


0.25 X 10"^ 


Nr 


1.023 


0.954 


0.834 


0.732 


Centrality [1.2 < \y\ < 2.2] 


0-20% 


20-40% 


40-60% 


60-92% 




0.91 X 10-3 


0.42 X 10-3 


0.19 X 10-3 


0.22 X 10-^ 


Nr 


1.023 


0.954 


0.834 


0.732 



ranging from — 20% to 60 — 92% are plotted in Fig. (4). The solid lines in the figures show 
the model-based results. 

Similarly, for Au + Au collisions at y'J/viv=200 GeV at RHIC, the values of hjj^, and Nr of 
equation (4) for different centrality regions are shown in Table 2 and the model-based results 
are plotted in Fig. 5. For calculating the values of Nr, we have used the values of < Npart > 
from [31]. The experimental results for the invariant yields of J/^ production as a function of 
transverse momenta at different centrality values ranging from — 20% to 60 — 92% and for the 
rapidities \y\ < 0.35 and 1.2 < \y\ < 2.2 respectively are taken from the PHENIX Collaboration 

m- 



3.3 Rapidity Distributions for Different reactions at y/s^N = ^OOGeV 

For the calculation of the rapidity distribution from the set of equations (1), (2), (3) and (4) we 

can make use of a standard relation as given below: 

dN _ r 1 <fN ^^^^ 
dy J 2ttpt dpxdy 

The rapidity distributions for the J/^'-production has now been reduced to a simple relation 
stated hereunder 

— — = ai exp(-0.23smh?/cm)- (12) 
dy 

The normalization factor oi depends on the centrality of the collisions and is obvious from the 
nature of the eqn.(l), eqn. (2), eqn.(3) eqn.(4) and eqn. (11). 

For p + p and d + Au collisions, the calculated rapidity distribution equations are 

^Ip+p^j/vp+x = 1-215 X 10-6exp(-0.23sinhyen^), (13) 
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Table 3: Values of ai x 10 for different centrality regions for Au + Au collisions 



Centrality 


0-20% 


20-40% 


40-60% 60-92% 


ai 


0.472 


0.033 


0.011 0.002 



and 

^U+A«^j/M = 7.025 X 10~^exp(-0.23sinhy,^), (14) 

In Fig. 6 and Fig. 7 we have plotted the rapidity distributions for J/^'-production in p + p and 
d + Au collisions respectively. Data in those figures are taken from Ref. ^28j and the lines show 
the theoretical outputs. 

To calculate the rapidity distribution for Au + Au collisions we are taking into account the 
different centrality regions by making use of Table 2 in addition to Eqs. (4) and (11). The values 
of ai are different for different centrality regions for Au + Au collisions and they are given in 
the Table 3. 

The solid lines in the Fig. 8 depict the theoretical plots of dN/dy vs. y for different centralities 
while the data for Au + Au collisions are taken from Ref. |30j . 

3.4 < pt > and < p^ > Values 

Next we attempt at deriving model-based expressions for both < px > and < p^ >. Of these 
twin observables the < pT > would be used by us in obtaining the ratios Raa-s for different 
participant numbers, Npart- Though, very strangely, there is so far no data on < px >"^/* even 
for p + p collisions at ISR to RHIC, we have chosen here to calculate and plot the nature of 

< px >'^/* vs. Npart for some collisions under study here for the following reasons: (i) in our 
opinion, < px >'^/* is more fundamental observable than < px (ii) for cases of almost 
all the secondaries produced in high energy collisions it is found that < px >^< Px > or 

< Px >'^^< Px > (iii) in our work that J/'^ behaves quite similarly as many other 
hadronic secondaries in nearly all the high energy interactions. So, plots oi < px > vs. Npart 
for Cu + Cu and Au + Au collisions are shown in Fig. 9 mostly in a predictive vein with the 
hope that it would be measured in future. And the excessive importance given to < pj^ >-^/* 
here by both the theorists and the experimentalists is just contextual. But, from physical 
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considerations, one must admit, < pT > is much more important than < pj. >. Furthermore, 
we also uphold the view that the theorists should not be guided by just the situational constructs 
of the experimentalists. 

The definition for average transverse momentum < px >"^/* is given below. 

npT (max) pd^a^ j„2 

(15) 

The theoretical plot of the average transverse momentum < px in Cu + Cu and Au + Au 

collisions at RHIC for different number of participating nucleons is shown in the Fig. 9. The 
theoretical calculations are done on the basis of uses of eqn.(4), Table 1 and Table 2. 

The QCD-oriented models believe that the transverse momentum distributions of the pro- 
duced J/^s carry valuable information about the mechanism for their absorption or disappear- 
ance manifested in the clearly falling nature of the production cross-section. The ideas based 
on the Standard Model (SM) predict that there should be a decrease of < > with centrality 
for sufficiently central collisions which is viewed as a consequence of the colour deconfinement 
'turnover effect '[32]. [33] . The present model-dependent behaviour of the observables < Px > 
as is defined by the undernoted relation (eqn. 16 given below). The [< >] is measured by 
RHIC and other experiments which is supposed to have a close relationship with what-is-called- 
to-be the " J/4' suppression", according to the QCD-points of view and to some QCD-oriented 
models. 

The expression for the average of the squared transverse momenta for any secondary (Q), by 
definition, is 

< p| >Q= ^ ^ (16) 

rPTimax) ^d^a'Jj 2 
JpT{min) dp-^ "T 

Using the above definition and putting the form of inclusive cross-section given by eqn. (4) into 
use, we calculate the < Pj^ > for the J/4' production in d + Au, Cu + Cu and Au + Au collisions 
collisions at RHIC energy ^snn = 200 GeV . The values of pTimin) and pximax) are taken 
here and 5 respectively [28], [30], [31] . 

We are interested not only in the absolute values of < > alone, but also in the centrality- 
dependence of them. Because, the conventional view holds the idea that < p^ > has a specific 
form of centrality-dependence, for which we present the graphical representations of Npart de- 
pendence of < p^ > for the rapidities \y\ < 0.35 and 1.2 < \y\ < 2.2 and for d + Au, Cu + Cu 
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Table 4: Comparison of < > values:: SCM vs. Hydro+J/^ [36] for different centrality bins 
and in mid-rapidity region in p+p and Au+Au collisions against the background of experimental 
values. 



Reactions Centrality Npart < p^ > {[GeV/c]^) < p^ > {[GeV/c]^) < p^ > {[GeV/c]^) 

(Expt.) (SCM) (Hydro+J/^) 



p + p 






2 


4.1 ±0.2 ±0.1 


3.65 


3.87 


Au + Au 


0- 


-20% 


280 


3.6 ±0.6 ±0.1 


3.64 


3.76 




20 


-40% 


140 


4.6 ±0.5 ±0.1 


3.66 


3.80 




40 


-60% 


60 


4.5 ±0.7 ±0.2 


3.55 


3.81 




60 


-92% 


14 


3.6 ±0.9 ±0.2 


3.43 


3.81 



and Au + Au collisions in Figures 10 (a) and (b) respectively. The experimental values are taken 
respectively from [28], [30], [3T] . 

Moreover, in Table 4, we present < p'^ > calculated from the two theoretical models, the SCM 
and the 'Hydro± J/'f' [36], with the experimental values [30] for p + p and Au + Au collisions 
at RHIC energies. 

3.5 The Nuclear Modifiacation Factor Raa 

There is yet another very important observable called nuclear modification factor (NMF), de- 
noted here by Raa which for the production of J/^ is defined by [30] 



<N,ouib)>d^N^j^^/dpTdy- 

the SCM-based results on NMFs for Cu + Cu and Au + Au collisions are deduced on the basis of 
Eqn.(4), Eqn.(7), Eqn.(8) and Table 1 , Table 2 and they are given by the undernoted relations 

^AA I c«+C«-J/M = 0.42p0.35 |y|<0.35, (18) 

Raa\cu+Cu^j/^+x = 0.38p?^35^ ^ 2 < \y\ < 2.2. (19) 

and 

Raa\au+Au^j/^+x = 0.36p!^i^ for \y\ < 0.35, (20) 
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Raa\au+Au-,j/^^+x = 0.26p?^l^ for 1.2 < \y\ < 2.2. (21) 

wherein the value of < Ncoii{b) > to be used is ~ 170.5 ib 11 [34J for Cu + Cu colhsions and 
for Au + Au colhsions it is taken as ~ 955.4 it 93.6 [35J. In Fig. 11, we plot Raa vs. pt for 
— 20% central region in Cu + Cu and Au + Au collisions. The solid lines in the figure show 
the SCM-based results against the experimental results [30], [31]. 

Next, for the analysis of centrality dependence as represented by (Npart), of the nuclear 
modification factor, (Raa) ^or Cu + Cu and Au + Au collisions at RHIC, we proceed in the 
direction with the help of eqn.(3), (4) and the calculated < pT >s [Fig. 5] for forward rapidity 
and different centralities. For Cu + Cu collisions, the data in Figure 12 are taken from [3D]. The 
solid line in that figure gives the SCM-based results. The dotted lines in Figure 12 show the 
hydrodynamical ('Hydro-|- J/^'') model-based calculations with and without nuclear absorption 
('hydro-|-j/Psi 1' and 'hydro-|-J/Psi 2' respectively) [37J. And for the Au + Au collisions, the 
data in Figures 13(a) and (b) are taken from |31j.|38j. The solid line in that figure gives the 
SCM-based results against the background of some other predictions or calculations [38], [39] 

3.6 LHC Prediction 

On the basis of the model-based calculations shown above, we now make the prediction for the 
invariant yields of J/^-mesons in the most central (0-10%) Ph + P6-collisions at the ^JsNN = 
5500 GeV i.e. at Large Hadron Collier (LHC) energy. The equation of the invariant yields of 
J/^* in the reaction of the type Ph + Pb J/'if + X fov the forward rapidity region would be 
of the form as given below 
1 (PN 

^5z/'^^|pb+Pb^j/*+x = 0.53 X 10-2p-i-i23gxp[-0.069(p^ + 9.61)] for 1.2 < \y\ < 2.2. 

(22) 

In Figure 14, we had presented the predicted nature of the yields for J/^ production in Pb + Pb 
collisions at LHC-energy by the solid curve. The dashed plot in the figure represents our RHIC 
calculations based on eqn.(4) and Table 2 in sect. 3.2. 

4 Summary and Conclusions 

In assessing the results achieved we must remember the limitations of the circumstances and 
the issues. It has to be appreciated that psions ( J/^'-mesons) are not as abundant as pions, the 
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commonest variety of the secondaries, the cross-sections for J/^'-production fall a few orders 
of magnitude compared to the pions. Besides, the uncertainties in the measurements of these 
psi-particles (resonance) are also much larger. So, in seeking the agreements of the model- 
based results with the measured data on J/^'-particles, one must be realistic. No model could 
probably claim a cent-percent agreement with data under such circumstances as pointed out 
above. For this reason, we have frequently compared here our model-based results with both 
the experimental data and the results obtained by one or two other competing models. There is 
yet another point very specific to J/^'-particles. The number of available and reliable data on 
psi-production in most cases is too limited. So, we will have to cut down our expectations about 
the degree of convergence between the calculations and the data and to arrive at the justified 
comments on any model-based work on J/^*. On this basis, we point out the following set of 
observations on the present work, (i) The agreement between our model-based (SCM-based) 
calculations and the pT-dependence of the invariant cross-sections on production J/^'-particles is 
modestly satisfactory, for mid and forward rapidities, moderate-p-r and central nuclear collisions, 
though there are some degrees of discrepancies for the other regions. This statement is valid on 
an overall basis for the figures from Fig. 2 - Fig. 5. (ii) The plots on studies of rapidity dependence 
are somewhat constrained by the very sparse data sets, for which we could make no strong or 
emphatic comment about the nature of agreements or disagreements. And this is a comment 
valid for the figures presented in Fig. 6 to Fig. 8. (iii) We discussed in the text given above 
somewhat extensively on the implications of the < p\ y^/^ and < px >'^/"^-plots; our results, 
based on the Figs. (9) and (10) and indicated by the Table-4 show fairly good agreement for 
< Px >"^/*-data obtained by nuclear collisions and also with 'Hydro+J/^' model [36j. (iv) 
And the results shown by Fig. (11)- Fig. (13) present a comprehensive study on the measured 
observable Raa and the model-based results, including the SCM-based plots in solid lines, (v) 
The Fig. 14 depicts our SCM-based predictions for the behaviour of the observable at LHC 
energy and in the LHC experiments at CERN. 

Thus, finally, we end up with a few conclusive remarks about the physical significance of 
the work: (i) As the standard model does not appear to be too promising for J/^I'-production 
behaviour, we have reasonably attempted a new one, surely with a certain degree of tentativeness. 
(ii) Neither the measurements for Cu + Cu interactions, nor the findings for Au + Au reactions 
at RHIC energies provide any final affirmation or clear indication of the QGP Physics. This 
is what is established by the very recent studies of both Tanenbaum [l],and Bhalerao, Gavai 
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[2]. (iv) The underlying mechanism that has been used here is not based on any 'gluonisation'- 
concepts; so this approach remains value-neutral to such theoretical inductions, (v) Nor does 
this mechanism view the J/^' particles as the products of the heavy quarks, for which the ideas 
of 'heavy quarks' cannot be accommodated in this model; Similarly the phrases like 'quark- 
gluon plasma' (QGP), 'anomalous suppression' etc. find no place in our work, (vi) But, the 
new mechanism for J/^'-production suggested by us, here incorporates the features of 'normal' 
suppression by the inclusion of the physical term of 'partonic rearrangement factor' yielding the 
normal suppression effect arising out mainly only at large transverse momentum values belonging 
to 'hard'-pr reactions, (vii) And for the 'hard' sector of hadronic collisions this model takes, 
without any exception, into account this 'constituent (or partonic) rearrangement factor' for 
production of all the secondaries (light or heavy)in a consistent manner, though certainly not on 
a uniform rate, which is certainly secondary-specific in nature. Thus, the 'globality property' of 
particle production at high energy (and at large pt) is upheld for the production of even J/^- 
mesons and so is not to be treated on a different basis in our approach, (viii) The anomalous 
observations [40] on the centrality dependence of the rapidity distributions very recently taken 
serious note of are kept consciously out of the purview of the present work and would be dealt 
with in a future work. 
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Pj(GeV/c) 

Figure 2: Plot of the invariant cross-section 
for J/^* production in proton-proton collisions at 
V^NN ~ "iOQGeV as function of px- The data 
points are from [26j. The solid curves show the 
SCM-based results. 
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Figure 3: Plot of the invariant cross-section for 
J/^* production in d -I- Au collisions at y^jvAf ~ 
200GeV as function of pr- The data points are 
from [28]. The solid curves show the SCM-based 
results. 
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Figure 4: Transverse momenta spectra at (a) \y\ < 0.35 and (b) 1.2 < \y\ < 2.2 for J/^* production 
in Cu + Cu collisions at Vs^^v ~ 2006*6^. The data are taken from [30J. The solid curves depict the 
SCM-based results. 
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Figure 5: Invariant spectra as function of pT for J /"^ production in Au + Au collisions at -v/stvtv ~ 
2Q{)GeV for (a) \y\ < 0.35 and (b) 1.2 < \y\ < 2.2 . The data are taken from [31 . The solid lines show 
the SCM-based results. 
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Figure 6: Plot of the rapidity distribution for J/^t production in proton-proton collisions at \/s^jv ~ 
200GeV as function of y. The data points are from 28 . The solid curves show the SCM-based results. 
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Fi gure 7: Rapidity distribution for J/^* production in d -I- Au collisions at \/sAf7V 
points are from [55] ■ The solid curves show the SCM-based results. 
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Figure 8: Plot of dn/dy vs. ?/cm-values for Au + Au collisions at y/sj^j^ = 2QQGeV for the centrality 
widths a)0 - 20%, b)20 - 40%, c)40 - 60% and d) 60 - 92%. The data points are from The solid 
line shows the SCM-based results. 
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Figure 9: Npart vs. calculated results on < px >-'/* based on SCM for (a) Cu + Cu and for (b) Au + Au. 
The calculations are done on the basis of eqn.(4), Table 1 and Table 2. No experimental data on 
< pt >-^/* are so far available. Calculated results are presented here in a predictive vein. 
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Figure 10: Plots on < p^ > versus Npart for J/^ production in d + Au, Cu + Cu and Au + Au collision 
in (a) mid and (b) forward rapidities. The data points are taken from [28], [30] and [31] and the hollow 
circles in the Figure show the SCM-based results. 
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Figure 11: Plot of Raa vs. p^-values for — 20% centrality in forward and mid-rapidities. The data 
points for Cu + Cu collisions (Fig. 8a) are taken from Ref.[3D] while those for Au + Au collisions (Fig. 
8b) are taken from Ref. ^l] . The solid lines show the SCM-based results. 
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Figure 12: Nature of iVport-dependence of the nuclear modification factor (NMF) denoted by Raa plots 
for Cu + Cu collisions at forward rapidity and at \/sj^j^ = 200CeV. The data are taken from Refer- 
ence [3D]. The SCM-based results are shown by the solid curved line in the Figure. The dashed curves 
show the predictions of the hydrodynamical model ('hydro-t- J/^') with (hydro+J/Psil)and without (hy- 
dro-|-J/Psi2) nuclear absorption [57] . 
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Figure 13: Nature of iVpart-dependence of the nuclear modification factor (NMF) denoted by Raa plots 
for Au + Au collisions at forward rapidity and at -\/sjvAr — SOOGel/. The data are taken from [3T|, [55] , 
The SCM-based result is compared with other theoretical results (a) Coalescence and Direct approaches 
[3H1 and (b)'hydro+J/^'' for Tj/^/Tc= 1.2 and 1.4(hydro+J/Psil, hydro+J/Psi2 respectively) [37], and 
PQCD from Ref. J9J in the Figures. 




Figure 14: SCM-based predicted plot on J/^f production in the most central Ph + Pb collisions at 
^/sWn ~ 5500 GeV i.e. at Large Hadron Collider (LHC) energy is shown by the solid curve. The dashed 
curve in the Figure represents our SCM-based plot for the most central collisions in Au + Au collisions 
at RHIC V^ivw = 200 GeV energy. 
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